****  itW#  Wrffeiqf 


ADAs/s# 


R-3698-AF 


Generalizations  of  Palm’s 
Theorem  and  Dyna-METRIC’s  j 
Demand  and  Pipeline 
Variability 


Manuel  J.  Carrillo 


January  1989 

A  Project  AIR  FORCE  report 
prepared  for  the 
United  States  Air  Force 


RAND 


i 

■, 

3 

1 


i 

> 

> 

3 

I 

i 

i 

i 


Approved  for  public  release;  distribution  unlimited 


04330-094^9 


best  available  copy 


'n  L*  C«'>rui)vnt.<o5rj- 


LaSSiFi  CATION  OF  This  »*GI  C<*h«n  Data  £n(w««) 


REPORT  DOCUMENTATION  PAGE 


5 


>  h- 


repo nr  number 

R-3698-AF 


title  fond  Subtitle) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  recipient’s  catalog  number 


3.  TYPE  OF  REPORT  A  PER1O0  COVERED 


Generalizations  of  Palm’s  Theorem  and  Dyna-METRIC  S  interim 
Demand  and  Pipeline  Variability 

PERFORMING  ORG.  REPORT  NUMBER 


AU  THORl'aJ 

Manuel  J.  Carrillo 


•.  contract  or  grant  numberTTJ 

F49620-86-C-0008 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

The  RAND  Corporation 
1700  Main  Street 
Santa  Monica,  CA  904Q6 

II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Requirements,  Programs  &  Studies  Group  (AF/RDQM) 
Ofc,  DCS/R&D  and  Acquisition 
Washington,  D.  C.  20030 


to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


12.  REPORT  OATS 

January  1989 

II.  NUMBER  OF  PAGES 

19 


1  MO»HTO«INGAClNCYNAM€*AOOWtSV7rdi7l5r«iilSiiirC«iMr«illit5o75«ij  IS.  SCCUAITY  CLASS,  (ol  thi*  export) 

Unclassified 

ts«.  oeClassification/oowngraoinc 
SCHEDULE 

Ui«VniviiTiwN  ffi.4  .‘A.4:  *®**—*=^ 

Approved  for  Public  Release;  Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (at  lha  akatraet  mi  farad  In  Blank  39.  It  dtttmranl  from  Report) 

No  Restrictions 


is.  Supplementary  notes 


'9.  KEY  RQRQS  (Contlmia  an  rararaa  si  da  il  nacaaaary  and  Idanilly  Or  block  number) 

Logistics  '  Poisson  Equation 

Inventory  Repair 

Models  Queueing  Theory  . 

Poisson  Density  Functions 
Distribution  Functions 

ao.  ABSTRACT  (Continue  on  smear aa  side  it  necaaaery  and  Idanilly  by  alack  number) 

J 

see  reverse  side 


oo ,  :«rn  1473 


SECURITY  CLASSIFICATION  of  THIS  page  .-Whan  Da, a  Entaead I 


Palm' s  theorem  is  a  useful  tool  in  modeling 
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nonhomogeneous  Poisson  input .  The  author 
shows  that  comparable  results  hold  for 
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PREFACE 


This  report  states  and  proves  the  theoretical  results  that  allow  the 
Dyna-METRIC  model  to  represent  dynamic  environments  and  parts 
failure  processes  with  a  range  of  variabilities.  These  results  are  exten¬ 
sions  to  the  field  of  queueing  theory  and  are  therefore  given  in  that 
context.  References  are  provided  for  the  reader  interested  in  applica¬ 
tions,  such  as  in  logistics  systems  of  the  military  services.  The  work 
presented  here  should  be  of  interest  to  logistics  modelers  as  well  as 
researchers  in  the  field  of  queueing  theory. 

The  research  leading  to  these  results  was  carried  out  in  the 
Resources  Management  Program  under  Project  AIR  FORCE. 


SUMMARY 


Palm’s  Theorem  is  first  reviewed  for  nonhomogeneous  Poisson 
arrivals  and  arbitrary  service  distributions  and  is  then  extended  to 
compound  Poisson  arrivals  and  arbitrary  initial  conditions.  For  these 
cases  it  can  represent  processes  where  either  the  number  of  demands  or 
the  repair  pipeline  evidences  extra-Poisson  variability.  For  many  ser¬ 
vice  distributions  of  interest,  the  resulting  probability  distributions  of 
the  number  in  the  queueing  system  have  algebraic  closed  form  solu¬ 
tions.  Therefore,  the  results  are  attractive  for  applications  in  modeling 
dynamic  environments  or  parts  failure  processes  having  a  wide  range  of 
variabilities.  One  such  application  is  in  the  field  of  logistics,  where 
these  results  play  a  major  role  in  the  Dyna-METRIC  model. 
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I.  INTRODUCTION 


Palm’s  Theorem  (Feller,  1968;  Palm,  1943)  for  the 
queue1  has  served  as  a  useful  tool  in  modeling  inventory  problems  in 
logistics  models  such  as  METRIC  (Sherbrooke,  1968)  and  Mod- 
METRIC  (Muckstadt,  1973).  However,  to  fit  its  limited  domain  of 
applicability,  time-dependent  customer  arrival  rates  have  been  approxi¬ 
mated  by  the  required  constant  rates.  And  even  when  the  time- 
dependent  behavior  is  modeled,  arrivals  with  large  variability  (as  mea¬ 
sured  by  variance-to-mean  ratios)  have  been  modeled  by  Poisson 
processes  (thus  having  a  variance-to-mean  ratio  of  unity).  Many  times 
this  has  caused  loss  of  accuracy  in  representing  the  problem  at  hand. 
Such  is  the  case  in  certain  military  logistics  applications  where  models 
need  to  represent  dynamic  environments  and  extra-Poisson  variation 
in  both  the  parts  failure  and  service  processes. 

This  report  reviews  the  extension  of  Palm’s  Theorem  for  time- 
dependent  arrival  rates  and  service  distributions  under  Poisson  input 
(Hillestad  and  Carrillo,  1980),  and  it  provides  further  extensions  to 
compound  Poisson  input.  Processes  can  now  be  modeled  where  the 
number  of  either  cumulative  arrivals  (demands)  or  customers  in  service 
(the  repair  pipeline)  has  variance-to-mean  ratios  greater  than  or  equal 
to  unity.  In  addition,  this  report  treats  a  range  of  initial  conditions  of 
the  queueing  system,  so  that  it  can  start  empty  or  at  a  fixed  or  random 
level. 

RAND’s  Dyna-METRIC  model2  can  represent  dynamic  environ¬ 
ments  for  logistics  support  systems  of  parts  failure  processes  having  a 
range  of  variabilities  and  initial  conditions.  Dyna-METRIC  has  found 
widespread  use  in  policy  analysis  and  evaluation  both  at  RAND  and  in 
the  Air  Force.3 

Section  II  introduces  the  nonhomogeneous  Poisson  queue  with  infi¬ 
nite  servers,  leading  to  a  generalization  of  Palm’s  Theorem  to  nonho¬ 
mogeneous  Poisson  input.  Section  III  shows  that  comparable  results 
hold  in  the  case  of  compound  Poisson  input.  Of  special  interest  is  the 
case  where  the  “compounding  distribution”  is  in  the  logarithmic  family, 

*In  queueing  notation,  the  M  stands  for  the  customer  interarrival  times  having  the 
exponential  or  memoryless  distribution,  the  G  stands  for  arbitrary  service  distribution, 
and  the  oo  stands  for  the  infinite  number  of  servers. 

2 Dyna-METRIC  is  a  model  developed  at  RAND  and  used  by  the  Air  Force  to  support 
policy  analysis  and  evaluation  of  multi-echelon  logistics  systems. 

3 Other  applications  of  the  results  given  here  are  suggested  in  Crawford  (1981);  Hille¬ 
stad  (1982);  Hillestad  and  Carrillo  (1980). 


1 


2 


which  allows  the  treatment  of  arrival  processes  with  variance-to-mean 
ratios  greater  than  unity.  Section  IV  relates  the  results  of  the  previous 
sections  to  two-echelon  repair  systems.  Finally,  Secs.  V  and  VI  cover 
various  initial  conditions  of  the  queueing  system. 


II.  THE  NONHOMOGENEOUS  POISSON  QUEUE 


DEFINITION  1.  Let  {N(t),  t  »  0}  be  the  counting  process  of 
arrivals  (e.g.,  failures  of  parts  of  a  given  type);  the  queueing  system 
starts  empty — N( 0)  =  0.  This  process  has  X(t)  as  intensity  function, 
X(t)  >:  0,  and  m(t)  as  mean  value  function, 

m(t)  -  E[N(t )]  =  /  Ms)  da, 
o 

where,  in  general,  £[7]  denotes  the  expected  value  of  a  random  vari¬ 
able  Y. 

DEFINITION  2.  Corresponding  to  an  arrival  that  occurs  at  time  s, 
let  the  positive  service  (repair)  time  random  variable  Y  have  an  arbi¬ 
trary  nonhomogeneous  distribution  depending  on  s,  so  that  the  proba¬ 
bility  that  Y  <y  is  P[Y  <  y]  =  G(s,  s  +  y).  The  service  times  of 
any  sequence  of  arrivals  are  assumed  to  be  independent,  and  the  distri¬ 
bution  G  is  assumed  to  have  a  finite  mean.  When  G(s,  s  +  y)  is 
homogeneous  (stationary),  it  will  be  denoted  by  F(y). 

Each  arrival  results  in  one  “customer”  in  the  nonhomogeneous  Pois¬ 
son  queue,  or  one  or  more  customers  in  the  nonhomogeneous  com¬ 
pound  Poisson  case  (Sec.  III).  This  report  characterizes  the  transient 
behavior  of  the  number  of  customers,  Z(t),  in  the  queueing  system, 
which  in  the  present  case  of  infinite  servers  is  also  the  number  in  ser¬ 
vice. 

As  a  point  of  departure,  a  generalization  of  Palm’s  Theorem  is  given 
in  a  form  that  describes  the  transient  behavior  of  the  M /G /o°  queue 
(see  Feller,  1968,  p.  460;  or  Palm,  1943). 

THEOREM  1  (Takacs,  1962,  p.  160): 

Suppose  arrivals  experience  independent,  identically  distributed  service 
times  from  a  stationary  distribution  F  as  described  in  Definitions  1 
and  2,  and  that  the  arrival  times  and  service  times  are  independent. 
For  the  special  case  of  constant  arrival  intensity  A(t )  -  c, 
c  >  0  for  t  >  0,  and  X(t)  -  0  otherwise  (homogeneous  Poisson  input 
starting  at  time  0),  the  number  of  arrivals  in  service,  Z(t),  has  a  Pois¬ 
son  distribution  with  mean 

A(t)  -  Jc[l  -  F(s )]  da. 
o 
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COROLLARY  1  (Takacs,  1962,  p.  160— Palm’s  Theorem): 

Under  the  assumptions  of  Theorem  1,  as  f  — *■  the  limiting  (steady 
state)  distribution  of  the  number  of  arrivals  in  service  is  Poisson  with 
mean 


lim  A(£)  =  cE[Y]. 

£— *  <*> 

Corollary  1  is  the  classical  form  of  Palm’s  Theorem.  Theorem  1 
generalizes  it  to  a  case  where  the  repair  distribution  is  homogeneous 
and  the  arrival  density  takes  a  jump  from  0  to  c  at  time  t0  -  0.  For 
processes  that  can  be  assumed  to  be  Poisson,  Corollary  1  has  provided 
the  basis  for  the  computation  of  spare  parts  levels  to  cover  the  pipeline 
of  parts  in  service. 

Because  a  sum  of  Poisson  random  variables  is  Poisson  with  mean 
equal  to  the  sum  of  the  means,  Theorem  1  suggests  that  homogeneity 
is  rather  unimportant— a  lot  of  nonhomogeneous  arrival  processes  can 
be  approximated  by  sums  of  several  arrival  processes  each  having  one 
jump.  In  fact  (as  will  be  shown  in  Theorem  2),  homogeneity  of  neither 
the  arrival  process  nor  the  repair  process  is  that  important. 

The  following  lemma  prepares  the  way  for  generalizing  Palm’s 
Theorem  to  nonhomogeneous  Poisson  arrivals. 

LEMMA  1  (Problem  13  in  Ross,  1970,  p.  29): 

For  the  arrivals  referred  to  in  Definition  1,  if  N(t)  -  n,  then  their  n 
arrival  times  have  the  same  distribution  as  the  order  statistics 
corresponding  to  n  independent  random  variables  each  having  the  dis¬ 
tribution 

H(x)  -  \m(x)/m(t)  for  0  <  x  <  t, 

[l  for  x  >  t. 

The  proof  is  along  the  lines  of  a  similar  theorem  for  the  stationary 
case  (as  given  in  Ross,  1970,  p.  17,  for  example).  It  is  now  possible  to 
prove  the  generalized  Palm’s  Theorem  for  nonhomogeneous  Poisson 
input. 

THEOREM  2  (Hillestad  and  Carrillo,  1980): 

Assume  that  the  arrivals  described  in  Definition  1  experience  indepen¬ 
dent  service  times  from  a  nonstationary  distribution  G  as  described  in 
Definition  2.  Furthermore,  assume  that  the  arrival  and  service  times 
are  independent.  Then,  the  number  of  arrivals  undergoing  service, 
Z(t),  has  a  Poisson  distribution  with  mean 
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A(t)  -  f  [1  -  G(s,  £)]X(s)  ds. 
o 

Hillestad  and  Carrillo  (1980)  provided  the  proof  of  this  theorem,  but 
it  is  also  given  here  for  completeness.  An  alternative  proof  appears  in 
Crawford  (1981).  For  a  historical  account  on  the  initial  conjecture  see 
Crawford  (1981);  Hillestad  (1982);  Hillestad  and  Carrillo  (1980). 

Proof: 

If  there  have  been  n  arrivals  by  time  t,  then  the  probability  that  k 
customers  are  still  in  service  at  time  t  is  given  by 

P[Z(t)  -  k  \N(t)  =  n]  -  (  g)p*(l  -P)"-*  for  k  -  0,  1,...,  n, 

zero  otherwise,  and  where 

p  -  [mlOl'1  J  [1  -  G(s,  £)]X(s)  ds 
o 

is  the  probability  that  an  arbitrary  arrival  is  still  in  service  (by  Lemma 
1).  Now,  conditioning  on  the  number  of  cumulative  arrivals  by  time  t 
yields 

P[Z(t)  -  *1  -  2  (  g  )  pHl  -  p)n~ke~m(t)mn(t)/ n! 


-  e-"*<‘>[ pm(t))k/k\  T  IU  -  p)m(£)]"-fe/(n  -  k)\ 

n-fe 


-  e-m<‘>A*(£)/fc!e(1-P>m('> 


-  e-A<‘>A k(t)/k\, 
using  the  fact  that  A(£)  -  pm{t). 

Because  Poisson-distributed  random  variables  have  the  same  mean 
and  variance,  the  variance-to-mean  ratios  of  both  the  number  of 
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arrivals,  X{t),  and  the  number  undergoing  service,  Z(t),  are  equal  to 
one. 

As  an  example,  consider  a  constant  intensity  function  A(£)  -  c, 
c  >  0,  and  exponential  service  times  with  mean  E[Y]  -  1/ft; 
Gfs,  8  +  y)  -  Ffy )  -  1  -  e-^.  Then  Aft)  -  (1  -  e~lU)c/fi.  As  a 
variation  of  this  example,  consider  the  case  where  the  servers  are  late 
and  become  operational  at  a  fixed  time  r  >  0.  This  leads  to  the  non¬ 
stationary  service  distribution 


Gfs,  t) 


0 

X  _  e-n(t-r) 
1  _ 


for  OssstiT, 
for  0  <  s  <  r  <  t, 
for  r  <  s  <  t. 


For  this  case, 


A  ft) 


ct 

cre-M-^  +  (1 


for  0  <  (  <  r, 
e-rtt-r))c  for  t  >  r  . 


Closed-form  algebraic  representations  of  Aft)  for  various  situations 
of  interest  are  given  in  Hillestad  (1982);  Hillestad  and  Carrillo  (1980). 
Their  ease  of  computation  make  them  attractive  for  applications. 


in.  THE  NONHOMOGENEOUS  POISSON  QUEUE 
WITH  COMPOUND  POISSON  INPUT 


The  results  presented  above  can  be  further  generalized  for  nonhomo- 
geneous  compound  Poisson  input.  The  last  section  used  the  term 
arrival  to  correspond  to  the  occurrence  of  a  part’s  failure.  The  more 
general  case  is  where  the  nth  arrival  results  in  the  failure  of  a  random 
number  Wn  of  parts.  In  queueing  parlance,  Wn  is  the  number  of  “cus¬ 
tomers”  in  an  arriving  batch,  and  Z(t)  is  the  total  number  of  custo¬ 
mers  in  service.  It  will  be  seen  later  that  the  Z(t)  process  is  useful  to 
represent  failure  processes  that  have  a  variance-to-mean  ratio  greater 
than  unity,  a  characteristic  rather  common  in  Air  Force  aircraft  com¬ 
ponent  removal  data.  Some  definitions  are  required. 

DEFINITION  3.  Let  Wn  be  the  batch  size  of  customers  in  the  nth 
arrival.  It  is  assumed  that  Wn,  n  =  0,  1,...  are  independent,  identically 
distributed  random  variables  having  the  common  compounding  distri¬ 
bution  {fj  -  P[W  -  j],j  -  0,  1,...}  having  mean  E[W\.  Further  it  is 
assumed  that  the  family  {Wn }  is  independent  of  the  arrival  process. 

DEFINITION  4.  Let  {X(t)>  t  >  0}  be  the  resulting  counting  pro¬ 
cess  of  customers,  where 

mt) 

X{t)  -  ±Wn, 

n-1 

and  JV(t)  is  as  given  in  Definition  1. 

As  is  apparent  from  a  conditioning  argument  on  N(t),  under  Defini¬ 
tions  1,  3,  and  4,  X(t)  has  the  compound  Poisson  distribution 

P[X(t)  -  k]  -  f.  fln'>e~m(t)mn{t)/n\  for  k  -  0,  1,..., 

n-0 

where  {fk  (n)}  is  the  n-fold  convolution  of  the  compounding  distribu¬ 
tion  {fk}  with  itself.  The  mean  and  variance  of  X(t)  are  given  by 

£[X(t)]  -  E[£[X(t)  |  N{t)]]  -  m(t)Em 
Var[X(f)]  -  £[Var[X(t)|N(t)]]  +  Var[E[X(t)  \  N{t)]] 
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-  Var[W]m(t)  +  E2[W]m(t) 

-  m(t)E[W*]. 

Therefore,  the  variance-to-mean  ratio  of  X(t)  is  E[W2]/E[W],  which 
is  only  dependent  on  characteristics  of  the  compounding  distribution 
and  not  on  t  or  m(t). 

It  is  now  possible  to  state  and  prove  the  generalization  of  Palm’s 
Theorem  for  nonhomogeneous  compound  Poisson  input. 

THEOREM  3: 

Let  the  arrivals,  the  compounding  distribution,  and  the  resulting 
nonhomogeneous  compound  Poisson  process  {X(£)i  i  2  0}  be  given  by 
Definitions  1,  3,  and  4,  respectively.  Also  assume  that  all  the  custo¬ 
mers  in  an  arriving  batch  have  equal  service  times  from  distribution  G 
in  Definition  2.  (Note  that  this  is  a  special  case  of  the  M /G /<*> 
queue  with  batch  arrivals  and  service.)  Then  the  distribution  of  the 
total  number  of  customers  in  service  is  compound  Poisson  with 

P[Z(t)  -  k]  -  2  /*0)e-A(t)A'(t)/;!  for  k  -  0, 1 . 

J-o 

and  where  the  parameter  A(t)  is  defined  as  in  Theorem  2  (and  thus, 
A (t)  is  here  the  expected  number  of  batches  in  service). 

Proof: 

By  Theorem  2,  the  number  of  batches  in  service  at  time  t  has  a 
Poisson  distribution  with  mean  A(t).  Therefore,  the  number  of  custo¬ 
mers  in  service,  Z(t),  is  nonhomogeneous  compound  Poisson  with  com¬ 
pounding  distribution  {/;,  j  -  0,  1,...}. 

COROLLARY  2: 

Under  the  assumptions  of  Theorem  3,  the  resulting  variance-to-mean 
ratio  of  the  number  of  customers  in  service,  Z(t),  is  also  given  by 
E[W2]/E{W]. 

COROLLARY  3: 

Consider  the  special  case  of  Theorem  3  having  constant  arrival  inten¬ 
sity  \(t)  -  c,t  a  0,  c  >0,  and  stationary  service  distribution  F  given 
in  Definition  2.  As  t  — * ►  <*>,  the  limiting  (steady  state)  distribution  of 
the  number  of  customers  in  service  is  compound  Poisson 
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Jim  P[Z(t)  -  k]  -  g  fkU)e-cElY](cE[Y]y/jl 
j-0 

The  result  given  in  Corollary  3  is  found  in  Feeney  and  Sherbrooke 
(1966). 

A  special  case  of  considerable  importance  in  applications  is  the  one 
in  which  the  compounding  distribution  of  the  input  is  the  logarithmic 
distribution. 

DEFINITION  5.  Let  W  have  the  logarithmic  probability  function1 

P[W  -  j] - (1  -  pV/0‘  logp),  0  <p  <  l;j  -  1,  2,... 

In  this  case,  E[W]  -  -(1/p  -  l)/logp  and  E[W2]/E[W)  -  1/p, 
which  implies  that  the  variance-to-mean  ratio  of  either  X(t)  or  Z(t) 
exceeds  one  and  can  take  on  arbitrarily  large  values  for  small  values  of 
P- 

It  is  well-known  (Feller,  1968,  p.  291)  that,  in  the  case  of  the  loga¬ 
rithmic  distribution,  X(t)  has  a  negative  binomial  distribution  with 
parameters  r(t)  -  -m(t)/log  p  and p.  Specializing  the  previous 
results  to  this  case  yields  the  following  theorem. 

THEOREM  4: 

Under  the  assumptions  of  Theorem  3,  if  the  compounding  distribution 
is  logarithmic  as  in  Definition  5,  then 

(a)  the  cumulative  number  of  customer  arrivals  X(£)  has  the 
negative  binomial  distribution  with  parameters 
r(t)  -  -m(t)/log p  and p,  and 

(b)  the  number  of  customers  in  service,  Z(t),  has  the  negative 
binomial  distribution  with  parameters  R(t)  -  -A(f)/logp 
and  p;  that  is, 

P[Z(t)  -  k]  -  (R(l)  \k  ~  Op^Hl  -p)*  for  k-0,1 . 

Proof: 

It  suffices  to  prove  case  (a),  since  the  proof  for  case  (b)  is  similar. 

Define  Ax  (s)  to  be  the  generating  function  of  the  probability  distri¬ 
bution  of  some  nonnegative  discrete  random  variable  X: 

‘The  parameter  p  in  this  logarithmic  distribution  is  unrelated  to  the  p  used  in  the 
proof  of  Theorem  2;  also,  log  p  is  the  natural  logarithm  of  p. 
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Axis)  -  g  siP[X  -  j ] 
j-  o 

For  the  nonhomogeneous  compound  Poisson  process  (see  Ross,  1970, 
p.  23,  for  example): 

Ax(t)($)  -  AN(t)(Aw(s))  -  exp[m(£)(-Aw(s)  -  1)] 


Therefore, 

4x<t)(s)  -  [e~l  exp  {log(l  -  (1  -  p)s)  /logp}]m(t), 
which  after  some  manipulation  leads  to 

A*U)(s)  -  [p/(l  -  a  -  p)s)]-'n(t)/l0*P, 


the  generating  function  of  the  required  negative  binomial  distribution. 

Section  I  referred  to  the  loss  of  accuracy  resulting  from  using  a 
model  with  Poisson  arrival  variability  (variance-to-mean  ratio  equal  to 
unity)  when  in  fact  the  arrivals  being  modeled  exhibited  greater  varia¬ 
bilities.  With  Theorem  4,  it  is  now  possible  to  see  the  effect  in  the 
number  of  customers  in  service  for  problems  with  compound  Poisson 
input  and  a  logarithmic  compounding  distribution.  It  suffices  to  com¬ 
pare  the  following  two  processes.  The  first,  with  variance-to-mean 
ratio  equal  to  unity,  is  a  nonhomogeneous  Poisson  process 
[N(t),  tiO)  with  mean  value  function  m(t),  as  in  Definition  1.  The 
second  is  a  nonhomogeneous  negative  binomial  process  with  variance- 
to-mean  ratio  1/p  >  1,  and  for  which  arrivals  have  mean  value  func¬ 
tion  -m(t){Iog  p)/(l/p  -  1)  so  that  both  processes  have  m(t)  as 
the  expected  number  of  cumulative  customers  arriving  by  time  t.  For 
the  first  process,  the  number  in  service  is  Poisson  with  mean  A (t), 
whereas  for  the  second,  the  number  in  service  is  negative  binomial  with 
parameters  A(t)/(l/p  -  1)  andp.  Statistics  for  these  two  distribu¬ 
tions  can  then  be  easily  obtained  and  compared  for  an  assessment  of 
the  magnitude  of  the  error  incurred  in  choosing  the  wrong  model. 


IV.  TWO-ECHELON  REPAIR  SYSTEM 
APPLICATIONS 


As  noted  at  the  beginning  of  this  report,  the  results  presented  here 
have  found  application  in  the  Dyna-METRIC  model.  Dyna-METRIC 
represents  the  failure  and  repair  of  components  resulting  from  the 
operation  of,  say,  a  fleet  of  aircraft.  In  that  context,  repair  of  com¬ 
ponents  occurs  at  either  at  the  base  level,  where  the  aircraft  operate,  or 
at  a  remote  service  depot;  this  is  known  as  a  two-echelon  repair  sys¬ 
tem. 

LEMMA  2  : 

Consider  the  nonhomogeneous  Poisson  process  given  in  Definition  1. 
Let  the  failed  part  corresponding  to  an  arrival  be  reparable  at  base 
level  with  fixed  probability  u  and  at  depot  level  with  probability 
(1  -  u),  independent  of  where  other  previous  failed  parts  are  repaired. 
Then  the  one  nonhomogeneous  Poisson  arrival  stream  with  mean  value 
function  m(t)  becomes  two  independent  nonhomogeneous  Poisson 
streams  with  mean  demand  functions  um(t)  for  the  base  and 
(1  -  u)m(t)  for  the  depot. 

The  proof  of  Lemma  2  is  along  the  lines  of  a  similar  result  for  the 
homogeneous  Poisson  case  (see,  for  example,  Ross,  1972,  p.  123). 

The  independence  of  the  base  and  depot  demand  streams  implies 
that  the  number  of  arrivals  undergoing  service  at  base  and  at  the  depot 
are  independent  and  that  the  distribution  of  their  sum  can  be  found  by 
a  convolution  of  their  distributions.  For  the  case  of  nonhomogeneous 
Poisson  input,  the  base  or  depot  number  of  arrivals  undergoing  service 
each  has  a  Poisson  distribution.  Therefore,  the  total  in  service  (base 
plus  depot)  has  a  Poisson  distribution  with  parameter  equal  to  the  sum 
of  the  parameters  of  the  base  and  depot  distributions. 

For  the  case  of  nonhomogeneous  compound  PoiBson  input,  the  base 
or  depot  number  of  customers  in  service  each  has  a  compound  Poisson 
distribution.  If  one  allows  the  base  to  have  a  compounding  distribution 
different  from  that  of  the  depot,  the  total  number  of  customers  in  ser¬ 
vice  is  not  necessarily  a  compound  Poisson  having  a  stationary  com¬ 
pounding  distribution. 

However,  when  base  and  depot  have  identical  logarithmic  com¬ 
pounding  distributions,  the  distribution  of  the  total  number  of 
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customers  in  service  is  compound  Poisson  with  the  same  compounding 
distribution— it  has  a  negative  binomial  distribution  with  parameters 
-  [A i(t)  +  A2(t)]/log  p  and  p,  where  the  subscripts  correspond  to  the 
base  and  depot  arrivals,  respectively. 

Sherbrooke  (1966)  dealt  with  two  homogeneous  independent  com¬ 
pound  Poisson  processes  (one  for  each  base  and  depot)  in  a  two- 
echelon  environment.  He  noted  that  if  the  two  compounding  distribu¬ 
tions  are  both  geometric  or  both  logarithmic,  then  the  distribution  of 
the  total  number  of  customers  in  repair  (base  plus  depot)  has  a  com¬ 
pounding  distribution  of  the  same  form  if  and  only  if  the  parameters  of 
both  compounding  distributions  are  identical.  Clearly,  that  result  has 
implications  for  the  case  of  nonhomogeneous  compound  Poisson  input. 
For  example,  if  the  compounding  distribution  is  logarithmic  and  either 
the  number  of  cumulative  customers  arrived  or  the  number  of  custo¬ 
mers  undergoing  service  at  the  base  had  a  variance-to-mean  ratio  dif¬ 
ferent  from  those  at  the  depot,  then  the  distribution  of  the  numbei  of 
customers  in  repair  at  base  and  depot  combined  would  no  longer  be 
negative  binomial.  Then  one  could  no  longer  take  advantage  of  the 
computational  simplicity  of  the  negative  binomial  distribution. 


V.  THE  NONHOMOGENEOUS  POISSON  QUEUE 
WITH  INITIAL  CONDITIONS 


This  section  relaxes  the  assumption  that  N  (0)  -  0 — that  is,  that  the 
queueing  system  starts  empty.  For  that  purpose,  consider  the  total 
number  of  arrivals  in  service  from  two  Poisson  processes,  one  whose 
arrivals  stop  at  t  -  0  and  another  that  starts  at  t  -  0.  The  results 
below  apply  to  the  nonhomogeneous  Poisson  queue,  but  similar  results 
for  the  nonhomogeneous  compound  Poisson  queue  are  obtained  in  Sec. 
VI.  To  arrive  at  steady  state  results  requires  restriction  to  the  case  of 
a  stationary  distribution  F  of  service  time  Y. 

DEFINITION  6.  Let  { JVi(t),  -T  <  t  <  »}  be  a  homogeneous  Pois¬ 
son  process  with  constant  intensity  function  Xi(t)  -  c  for  - 71  <  !  <0, 
and  zero  otherwise.  Furthermore,  this  system  starts  empty — 
N^-T)  «  0. 

DEFINITION  7.  Also  let  {JV2(t),  t  >  0}  be  a  nonhomogeneous 
Poisson  process  with  intensity  function  X2(t),  comparable  to  the  pro¬ 
cess  in  Definition  1  and  independent  of  the  process  in  Definition  6. 

The  arrivals  of  these  two  processes  have  service  times  from  a  sta¬ 
tionary  distribution  F  as  in  Definition  2,  consistent  with  the  idea  of 
having  a  common  set  of  an  infinite  number  of  servers.  Then 
Zi(t)  and  Z2(t),  the  number  of  arrivals  in  service  from  the  first  and 
second  processes,  are  independent.  It  follows  that  the  total  number  of 
arrivals  in  service,  Z(t)  -  Z\(t)  +  Z2(t),  has  a  distribution  given  by 
the  convolution  of  the  corresponding  distributions  of  the  two  processes. 

This  problem  arises  in  military  logistics  applications.  For  example, 
an  aircraft  fleet  may  be  flying  at  an  approximately  constant  rate  in 
peacetime,  but  then  a  military  condition  arises,  at  t  -  0,  where  the 
intensity  of  flying  increases  and  changes  over  time.  The  operation  of 
the  fleet  of  aircraft  causes  aircraft  parts  to  fail  and  require  repair;  in 
such  circumstances,  these  queueing  models  could  represent  the  number 
of  each  part  type  undergoing  repair  at  some  positive  time  t  (Hillestad 
and  Carrillo,  1980). 

Corollaries  4  and  5  below  pertain  to  two  types  of  initial  conditions 
for  one  nonhomogeneous  Poisson  process  that  is  preceded  by  another 
that  is  homogeneous  Poisson.  One  set  of  initial  conditions  assumes 
knowledge  of  the  number  of  arrivals  in  service  at  t  -  0,  whereas  the 
other  assumes  that  at  t  -  0  the  number  in  service  has  a  limiting 
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(steady  state)  distribution.  These  results  are  applicable  to  all  (measur¬ 
able)  stationary  service  distributions  with  a  finite  mean. 

LEMMA  3: 

Consider  the  arrival  process  in  Definition  6  with  stationary  service  dis¬ 
tribution  F  in  Definition  2.  It  is  further  assumed  that  the  arrival  and 
service  processes  are  independent.  If  Zi(0)  -  n,  then,  as  T  — ►  <»,  the 
limiting  distribution  of  the  number  of  arrivals  still  in  service  at  time 
t,  Z\{t),  t  >  0,  is  binomial  with  parameters  n  and  [1  -  F*(t)],  where 

F'(t)  -  [tffY]]-1  j  [1  -  F(s)]  ds 
o 

is  the  equilibnum  distribution  of  F  (see  Ross,  1970,  p.  47). 

Proof: 

The  proof  of  Lemma  3  is  based  on  a  result  found  in  Takacs  (1962, 
Theorem  2,  p.  161)  basically  stating  that,  given  Zi(0)  -  n  and  steady 
state,  the  probability  that  a  typical  customer  in  service  at  time  0  will 
depart  by  time  t  is  given  by  the  equilibrium  distribution  of  F. 

Note  that  the  equilibrium  distribution  for  an  exponential  is  the  same 
exponential.  As  a  second  example,  consider  the  degenerate  (deter¬ 
ministic)  distribution 

F(y)  -  Ofory  <  E[Y)  and  F(y)  =  lfory  s  E[Y]. 

The  equilibrium  distribution  in  this  case  is  uniform  with 


F*(y)  - 


0  for  y  <  0, 
y  for  0  <  y  <  E[Y], 
1  for  y  2;  E[Y]. 


Furthermore,  note  that  if  F  is  exponential,  Lemma  3  holds  for  any 
T  >  0:  the  steady  state  requirement  is  not  needed  because  of  the 
memoryless  property  of  the  exponential  distribution. 


COROLLARY  4: 

Consider  the  combined  processes  {Ni{t),  -T  <  t  <  <»}  and 
{7V2(£),  t>0}  given  in  Definitions  6  and  7,  along  with  stationary  ser¬ 
vice  distribution  F  as  in  Definition  2.  Again,  the  arrival  process  and 
the  service  times  are  assumed  independent.  If  Z(t)  -  Zx(t)  +  Z2(t), 
then,  as  T  — ►  oo,  the  limiting  distribution  of  P[Z(t )  -  k  |  Zi(0)  -  n]  is 
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the  convolution  of  the  binomial  distribution  in  Lemma  3  (for  Zi(t)) 
with  the  Poisson  distribution  from  Theorem  2  (for  Z2(t)  using  X2(t)). 

As  an  example,  consider  the  special  case  of  constant  input  rate 
X2(f)  -  c  and  exponential  service  distribution  F  with  mean 
i?[Y]  -  1/m  so  that  A2(t)  *  (l-e_*“)c/M-  Corollary  4  then  leads  to 
the  transition  probabilities  P,y(t)  for  the  M/M queue  as  given  in 
Saaty  (1961,  p.  100);  Feller  (1968,  p.  481): 

Pin(t)  -  lim  P[Z(t )  -  n  |  Z{ 0)  -  i]  * 

T— °° 

£(  -  e-^)*'*  e-Al(t)[A2(t)]n-V(n  -  k)\ 

fc-o' 

for  i  -  0,  1,...;  n  =  0,  1,...,  t;  t  >  0. 

For  the  case  n  >  i,  the  above  expression  is  valid  if  the  upper  limit  of 
the  summation  index  is  changed  to  i. 

Now  we  turn  to  different  initial  conditions,  from  one  of  known 
number  of  customers  in  the  system  at  t  -  0  to  an  assumption  of  ran¬ 
dom  traffic. 

LEMMA  4: 

Consider  (Nx(t),  -T  st  <»)  in  Definition  6  and  the  stationary  ser¬ 
vice  distribution  F  in  Definition  2.  Assume  that  the  arrival  process 
and  the  service  times  are  independent.  Then,  for  t  >  0, 

A i(t)  -  *  f  c[l  -  F(s )]  ds, 

t 

and,  as  T  — ►  oo,  the  limiting  distribution  of  the  number  of  arrivals  in 
service  Z\(t)  is  Poisson  with  mean 

lim  Ai(t)  -  c£;[Y){l  -  F’(t)} 

T- 00 

Proof: 

By  Theorem  1, 

0  t  +  T 

Ai(t)  -  J* c [1  —  F(t  -  s)J  ds  -  j  c[l  -  F(s )]  ds. 

-T  t 
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Therefore, 


t 


thus  completing  the  proof. 


It  is  of  interest  that  F'(t),  the  equilibrium  distribution  of  F,  plays  a 
role  in  the  limiting  distributions  of  both  types  of  initial  conditions 
treated  in  this  section.  Also,  for  t  -  0,  the  above  result  shows  that  the 
steady-state  distribution  of  the  number  of  arrivals  in  service  is  Poisson 
with  parameter  cl£[  Y],  in  agreement  with  Corollary  1. 

As  an  example,  if  F  is  exponential  with  mean  E[Y)  -  l/n,  -hen 

A i(t)  -  e (1  -  e~T/l,)c/n  and 


lim  A\(t)  -  e_"‘c/M. 

If  F  is  the  degenerate  (deterministic  time)  distribution, 


Ai(f)  - 


cT 

c(E[Y]  -  t) 
0 


for  0<t<  E[Y]  -  T, 

for  £[Y)  -  T  <t  £  E[Y], 

for  t  >  E[Y]; 


lim  Ax(t)  -  (  Y]  -t)  for  0<(s  E[Y], 
r-.  \  0  for  t  >  E[Y]. 


COROLLARY  5: 

Consider  the  assumptions  of  Corollary  4.  For  the  combined  processes 
{Ni  (f),  -  T  si  <»)  and  {N2(t),  t  a  0),  the  distribution  of  the 
number  of  arrivals  in  service,  Z(t),  is  Poisson.  For 
t  <  0,  7At)  -  Zj(t);  for  t  >  0,  Z(f)  -  Zi(t)  +  Z2(t).  In  either  case, 
the  corresponding  A(t)  are  found  by  appropriate  applications  of 
Theorem  2,  and  the  limiting  form  of  Ai(t)  as  T  —*  <»  is  given  by 
Lemma  4. 

Clearly,  Corollary  5  is  quite  attractive  for  modeling  because  of  its 
computational  simplicity. 


VI.  THE  NONHOMOGENEOUS  COMPOUND 
POISSON  QUEUE  WITH 
INITIAL  CONDITIONS 


The  results  of  Sec.  V  are  valid  for  the  nonhomogeneous  Poisson 
processes  that  govern  the  arrival  of  batches  in  nonhomogeneous  com¬ 
pound  Poisson  processes.  In  this  case,  the  total  number  of  customers 
in  service  (from  the  arriving  batches),  Z(t)  =  Z i(£)  +  Zz(t),  has  a  dis¬ 
tribution  given  by  the  convolution  of  the  corresponding  distributions  of 
the  two  processes. 

The  following  corresponds  to  Lemma  3. 

LEMMA  5: 

Let  the  Poisson  process  and  compounding  distribution  in  Definitions  6 
and  3,  respectively,  specify  a  compound  Poisson  process.  Also,  let 
resulting  customers  be  provided  service  from  a  stationary  service  distri¬ 
bution  F  as  in  Definition  2;  the  arrival  process  is  assumed  independent 
of  the  service  times.  Then,  for  t  >  0  and  as  T  —  a  limiting  distri¬ 
bution  for  the  number  of  customers  still  in  service  at  time  t  is  given  by 

lim  P[Zi(t)  -  k  \Zi(0)  -  n]  = 

T— - 


{lim  P(ZX( 0)  -  n]}"1 

r— oo 

2  2  /w(  -  hY->f"e-'W\(cE[Y)V/i\ 

i-0j-0 


for  h  -  1  -  F*(t);  k,  n  -  0,  1,... 

If  the  compounding  distribution  is  logarithmic,  then,  as  T  —  »,  the 
limiting  distribution  of  the  number  of  customers  in  service  at  t  -  0  is 
negative  binomial  with  parameters  -c£[Y]/ log  p  and  p. 

Proof: 

Let  B(t)  be  the  number  of  arrival  batches  in  service  at  time  t  >  0. 
Then  conditioning  on  B  (0)  -  i  yields 
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PlZ.it)  =  k  |Zi(0)  =  n]  = 

2  P[Z,(t)  =  k  |Z,(0)  =  n,B(0)  =  i]P[B(0)  =  i  |Zi(0)  =  n]. 

i-0 

But  because  of  Bayes’  formula, 
lim  P[B (0)  =  i  |Zi(0)  -  ra]  - 


{lim  P[Z!(0)  =  n)}-lf^e-cE^\cE[Y)y/i\ 

T- oc 

However, 

P [Z j(t )  =  |Zj(0)  =  n,  B(0)  -  i]  - 

2  P[Z!(t)  -  A,  fltt)  =  ;  !  Zj(0)  =  n,  B(0)  >  i), 

j-0 

which  leads  to  the  desired  result. 

The  following  corresponds  to  Corollary  5. 

COROLLARY  6: 

Consider  the  two  processes  {N.(t),  -T  si  <»)  and  {iV2(t),  t  >  0}  in 
Definitions  6  and  7  having  the  identical  logarithmic  compounding  dis¬ 
tribution  as  specified  in  Definitions  3  and  5.  The  stationary  service 
distribution  F  is  given  in  Definition  2,  and  the  service  process  is 
independent  of  the  arrivals.  Then,  for  the  combined  process,  the  dis¬ 
tribution  of  the  number  Z(t)  of  customers  in  service  (from  arriving 
batches)  is  negative  binomial.  For  t  <  0,  Z(t)  -  Zj(t);  for  t  a  0, 
Z(t)  =  Z .it)  +  Z2(t).  In  either  case,  the  corresponding  A(t)  are  found 
by  appropriate  applications  of  Theorem  4  or  Lemma  4.  Thus  Z(t), 
t  >  0,  has  a  negative  binomial  distribution  with  parameters 
—  [Al(t )  +  A2(t)]/log  p  andp. 
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